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LASER SYSTEM USING ULTRA-SHORT LASER PULSES 

CROSS REFERENCE TO RELATED APPLICATIONS 
[0001] This is a continuation-in-part of U.S. Patent application Serial 
No. 10/265,211, filed October 4, 2002, which is a continuation-in-part of PCT 
Serial No. PCT/US02/02548, filed on January 28, 2002, which claims priority to 
U.S. provisional application Serial No. 60/265,133, filed January 30, 2001, which 
are all incorporated by reference herein. 

BACKGROUND 

[0002] The present invention generally relates to a laser system and 
more particularly to a laser system using ultra short laser pulses and a pulse 
shaper. 

[0003] Commercially practical femtosecond lasers have been 
unavailable until recently. For example, lasers which can generate 10 
femtosecond or less laser pulse durations have traditionally been extremely 
expensive, required unrealistically high electrical energy consumption (for 
extensive cooling, by way of example) and depended on laser dyes that had to 
be replenished every month thereby leading to commercial impracticality. The 
efficiency of sub-10 femtosecond lasers was not practical until the year 2000 
because of the prior need for dyes and flash lamps instead of YAG and Ti: 
Sapphire crystals pumped by light or laser emitting diodes. 
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[0004] Ultrashort pulses are prone to suffer phase distortions as they 
propagate through or reflect from optics because of their broad bandwidth. 
There have been recent experimental attempts to shape the phase of ultrashort 
pulses since shaped pulses have been shown to increase the yield of certain 
chemical reactions and multiphoton excitation. Pulse shaping methods allowing 
synthesis of complex femtosecond optical waveforms are known. As usually 
practiced, the output waveform is determined by the Fourier transform (FT) of a 
spatial pattern transferred by a mask or a modulator array onto the dispersed 
optical spectrum. As an example, the FT optical pulse shaper is widely used for 
synthesis of complex femtosecond waveforms. In this geometry, the temporal 
profile of the output waveform is given by the FT of the mask pattern that is 
transferred onto the optical frequency spectrum of the pulse. FT pulse shaping 
was first demonstrated for use in simple pulses of tens of picoseconds in 
duration. Pulse shaping was then extended to the sub-100 femtosecond (fs) (10' 
15 second) time scale and demonstrated highly structured waveforms using 
microlithographically patterned pulse shaping masks. The introduction of liquid 
crystal modulator arrays and acousto-optic (A/O) modulators into FT pulse 
shapers led to computer programmable pulse shaping, with millisecond and 
microsecond reprogramming times, respectively, and widespread adoption of this 
technique. 

[0005] These shaped pulses require a very large data set and in many 
cases, complex learning calculations for determining the pulse shaping 
characteristics for a particular application. The optimal pulse for the particular 



2 



6550-000057/CPE 

application is not known in advance. Since the variation shape of the possible 
pulse shapes is huge, scanning the entire parameter space is impossible and as 
such the optimized pulse shape could not have been predicted by theory. For a 
pulse shaper with N pixels, one can generate (P *A) N shaped pulses, where P 
and A are the number of different phases and amplitudes a pixel can take. If it is 
assumed 100 pixels, each taking 10 different amplitude values and 100 different 
phase values, the number of different pulses is of order of magnitude 10 300 . This 
dataset is extremely large, therefore, while in principle, the field exists to achieve 
the desired photonic transformation or excitation, finding it is a great challenge. It 
would be desirable for a system to control ultrashort pulses with a smaller 
dataset and operable to generate very complex pulse shapes that are optimal for 
the particular application. 

SUMMARY OF THE INVENTION 
[0006] In accordance with the present invention, a laser system using 
ultrashort laser pulses is provided. In another aspect of the present invention, 
the system includes a laser, pulse shaper and detection device. A further aspect 
of the present invention employs a femtosecond laser and binary pulse shaping 
(BPS). Still another aspect of the present invention uses a laser beam pulse, a 
pulse shaper and a SHG crystal. In yet another aspect of the present invention, 
a multiphoton intrapulse interference phase scan (hereinafter "MIIPS") system 
and method characterize the spectral phase of femtosecond laser pulses. In 
another aspect of the present invention, a system employs electromagnetic pulse 
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shaping design to take advantage of multiphoton intrapulse interference. Fiber 
optic communication systems, photodynamic therapy, functional imaging, and 
pulse characterization tests use the laser system with additional aspects of the 
present invention. 

[0007] The present invention also improves the encoding-decoding 
functionality of pulses by adding considerably more information to each pulse by 
obtaining the entire phase function directly from a phase scan. Intrapulse 
interferences of the present invention causes self separation (for example, 
inherent communication signal routing address differentiation) thereby allowing 
use of inexpensive receivers in an asynchronous manner, in other words, without 
the need for synchronous detection such as by traditional autocorrelation or 
interferometers. The control of nonlinear optical processes using multiphoton 
intrapulse interference can be applied in diverse fields such as photochemistry, 
communications, and medicine. 

[0008] Further areas of applicability of the present invention will 

become apparent from the detailed description provided hereinafter. It should be 
understood that the detailed description and specific examples, while indicating 
the preferred embodiment of the invention, are intended for purposes of 
illustration only and are not intended to limit the scope of the invention. 

BRIEF DESCRIPTION OF THE DRAWINGS 
[0009] The present invention will become more fully understood 

from the detailed description and the accompanying drawings, wherein: 
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[0010] Figure 1 is a diagrammatic view showing a preferred 
embodiment of a laser system of the present invention; 

[0011] Figures 2a-2d illustrate phase scans from MIIPS as compared 
to FROG; 

[0012] Figure 3 is a graph showing phase scans; 

[0013] Figures 4a-4c are graphs showing phase scans; 

[0014] Figure 5 is a diagrammatic view showing a preferred 
embodiment of the present invention that employs MIIPS; 

[0015] Figure 6a is a flowchart of computer software that is used in 
some of the preferred embodiments of the invention; 

[0016] Figure 6b is a flow chart of computer software that is used in an 
automated spectral phase determination of some of the preferred embodiments 
of the present invention; 

[0017] Figure 7 is a diagrammatic view of a preferred embodiment of 
the present invention applied to multi-photon microscopy; 

[0018] Figure 8 is a graphical illustration of the contrast ratio C that can 
be obtained by two-photon excitation using BPS; 

[0019] Figure 9 is a graphical illustration of a pseudo-random number 
generated BPS using prime numbers; 

[0020] Figures 10a and 10b demonstrate experimental results and 
theoretical predictions of a method based on amplitude masking not using BPS; 
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[0021] Figures 1 1 a-1 1 b are graphical representations of the pseudo- 
random number generation BPS and an output after second harmonic 
generation; 

[0022] Figures 12a-12d a graphical illustration showing experimental 
results and theoretical results from BPS employed in the present invention, 
before and after optimization with an evolutionary learning ; 

[0023] Figure 13 is a diagrammatic view showing embodiments of the 
present invention applied to optical coherent tomography and photo dynamic 
therapy; 

[0024] Figure 14 is a diagrammatic view showing another preferred 
embodiment of the invention applied to spectroscopy; 

[0025] Figure 15 is a diagrammatic view showing a preferred 
embodiment system of the present invention using BPS applied to optical 
coherent tomography and photodynamic therapy; 

[0026] Figure 16 is a diagrammatic view showing an another preferred 
embodiment of the invention applied to communications; 

[0027] Figure 17 is a diagrammatic view showing another preferred 
embodiment of the present invention applied to communications; 

[0028] Figures 18a and 18b are diagrammatic views showing 
components in an alternative preferred embodiment applied to communications; 

[0029] Figures 19a and 19b are experimental results and theoretical 
predictions of selective two photon excitation; 
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[0030] Figures 20a-20f are experimental results of pH-sensitive two 
photon microscopy; 

[0031] Figures 21a-21d are experimental results of selective two 
photon microscopy sample doped with different fluorescent probes. 

[0032] Figures 22a-22d are experimental results of selective two 
photon microscopy of microspheres; and 

[0033] Figure 23 is a perspective view showing a preferred 
embodiment of a fixed, two dimensional shaper employed in the present 
invention. 

DETAILED DESCRIPTION OF THE PREFERRED EMBODIMENTS 
[0034] The following description of the preferred embodiments is 
merely exemplary in nature and is in no way intended to limit the invention, its 
application, or uses. 

[0035] The present invention provides methods and apparatus for a 
laser system using ultra short pulses. Pulse shaping essentially involves control 
over the amplitude, phase, frequency and/or inter-pulse separation. Complex 
pulse shaping aims to control one or more of the above-mentioned parameters in 
a programmable manner, such that the user has complete control. In other 
words, complex pulse shaping allows generation of complicated ultrafast optical 
waveforms according to user specification. 

[0036] Coherent control is the ability to control the dynamics at various 
stages of a process as it evolves under the effect of a coherent source. Many of 



7 



6550-000057/CPE 

the frequencies constituting the ultrafast pulse can simultaneously excite many 
coherent transitions to the excited states, and a capability to manipulate them 
with the shaped pulses lead to the results. The term "coherent control" refers to 
the manipulation of molecular states coherently and thus avoids the issue of 
uncertainty principle for situations involving ultrafast timescales taking advantage 
of constructive and destructive interference. 

[0037] One of the easiest available pulse modulation schemes is 
frequency chirping. "Chirping" essentially refers to the process of arranging the 
frequency components in a laser pulse with certain phase ordering. Linear 
"ordering" can be easily achieved by dispersing the ultrafast pulses through a 
pair of grating or through pulse propagation in optical fiber. Due to the uncertainty 
principle, this "ordering" of frequency components results in the lengthening of an 
otherwise bandwidth-limited ultrafast pulse. Typically, linear frequency chirping is 
the most often used shaped pulses, although higher-order chirp generation and 
some of their recent experimental uses have also been predicted. 

Multiphoton Intrapulse Interference 

[0038] A multiphoton intrapulse interference phase scan (hereinafter 
"MIIPS") system and method of the present invention characterize the spectral 
phase of femtosecond laser pulses. The phase across the spectrum of an 
ultrafast pulse can affect the multiphoton process in a number of ways. Phase 
can increase the pulse length and hence reduce the peak intensity of the pulse, 
thereby preventing saturation, a common result under high intensity excitation. 
Phase can also be used to synchronize changes in the electric field with 
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intramolecular wave packet dynamics. This idea has been explored theoretically 
and invoked to explain some gas and condensed phase experiments. Finally, 
phase can be used to cause interference in the way multiple frequencies 
combine to achieve multiphoton excitation. This process also known as 
multiphoton intrapulse interference (Mil) has been demonstrated under sharp 
resonance conditions such as transitions in isolated atoms (See Broers, B., et al., 
"Large interference effects of small chirp observed in 2-photon absorption," Opt. 
Commun. 1992, 91, 57; and Broers, B., et al. "Diffraction and focusing of spectral 
energy in multiphoton processes," Phys. flev. A 1992, 46, 2749.) 

[0039] The technique of MM and its application to control multiphoton 
processes is based on rationally designing an electric field required to achieve an 
articular target with a minimum number of parameters. The method is based on 
calculating the amplitude of the nth-order electric field and comparing it to the 
absorption spectrum of the molecules being controlled. This provides a strong 
physical understanding of the control process, which can be very useful in the 
interpretation of experiments where the field is optimized by computer programs 
based on evolutionary learning s or similar methods. 

[0040] Two-photon transitions can focus the energy from an ultrafast 
pulse into a narrow frequency distribution; just like Fresnel diffraction can be 
used to construct a focusing lens. Conceptually, Mil takes advantage of the 
interference term that is associated with the phase of each frequency v/' within 
the pulse that contributes to the multiphoton process and can enhance or what 
may be as valuable, suppress a multiphoton transition. The effective electric field 
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that drives the two-photon process through the induced (nonlinear) polarization is 
proportional to E 2 (f) (in the absence of intermediate resonance at the one-photon 
level). Its Fourier transform £< 2) (v) determines the frequency response at the two- 
photon level. 

[0041] A new method known as multiphoton intrapulse interference 
phase scan (Ml I PS) is capable of both pulse characterization and compensation 
of subsequent pulses. Within minutes, the pulses are characterized and 
compensated to yield transform-limited (TL) or user-specified shaped pulses at 
the sample. This capability is extremely practical and can be incorporated in any 
laser setup. 

[0042] MIIPS is a single-beam method that does not require an 
interferometer (see Figure 1). To make a precise and accurate measurement of 
the spectral phase using MIIPS, a known phase delay is imposed on the 
frequencies that make up the pulse using a calibrated pulse shaper. The pulse 
shaper essentially behaves as two back-to-back spectrometers. In one 
embodiment, the pulse is dispersed with a prism and collimated with a 200-mm 
cylindrical mirror. At the Fourier plane, where all the frequencies are isolated, 
their phases are manipulated by a computer-controlled LCD spatial light 
modulator (SLM). The SLM applies the reference phase function to the input 
pulse, and the resulting pulse is then reconstituted to the time domain by a 
second cylindrical mirror and prism. The SLM can be updated every pulse 
(presently limited to 1 kHz). The LCD has a 250-ms response time, so in 
principle it can be updated at 4 kHz. The output beam is analyzed by placing a 
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0.01-mm-thick beta barium borate crystal for second-harmonic generation (SHG) 
in its path, usually at the place where optimum pulses are required. The use of 
the second harmonic is critical to the method. In a sense, the pulse 
autocorrelates itself at the SHG crystal. For each reference phase function that is 
introduced by the computer-controlled SLM, the output spectrum from the SHG is 
dispersed in a spectrometer and recorded. 

[0043] Pulse characterization involves the introduction of a reference 
phase-modulation function of the form <p= cos (y Q - 5), where °c. is the 
magnitude of the phase delay, y is the periodicity Q is the frequency detuning 
from the carrier frequency of the pulse, and 5 is the position in the spectrum at 
which the cosine function is equal to one. The reference phase function, with 
typical values « = 2n, and y= pulse duration, is programmed into the SLM and 
scanned for different values of 5 ranging from 0 to 2 n. For each value of 5, the 
spectrum of the frequency- doubled pulse changes, achieving a maximum in the 
spectral region over which the SLM compensates for the phase distortions. The 
Ml I PS trace corresponds to the collection of spectra as a function of 5 (see 
Figure 2). MllPS-generated trace of wavelength as a function of 6 shows 
changes in the SHG spectrum of the laser pulse intensity. In general, the 
distance between the diagonal features is proportional to linear chirp and the 
angular deviation is proportional to quadratic chirp. Computer analysis of the 
trace is used to retrieve the spectral phase of the input pulse (a). The FROG 
(frequency resolved optical gating herein known as "FROG") trace clearly shows the 
spectral phase distortion of the pulse by its deviation from an oval (b). After three 
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iterations of characterization and compensation, the output pulses are transform- 
limited as evidenced by the parallel features in the MIIPS data (c) and the oval 
feature in the SHG-FROG (d). 

[0044] Qualitatively, the distance between the diagonal features 
determines linear chirp while the angle between the features determines the 
quadratic chirp. The full quantitative determination of the spectral phase by 
integration can be obtained. Once the MIIPS system has characterized the pulse 
and retrieved the phase distortions inherent to the pulses, it can use that 
information to drive the SLM such that it compensates for the distortions. The 
first step in compensation is to take the phase determined from the first scan and 
program it into the SLM with a negative sign so that it subtracts the distortions. 
The system carries out a new phase scan to determine the remaining spectral 
phase modulation (usually about 10% of the original). Typically, three such 
iterations will yield transform-limited pulses. Because the laser is not focused in 
the pulse shaper, the method can be used with pulses that are relatively high in 
energy. Pulses ranging from about 100 pJ to about 1 mJ and pulse durations 
from less than 5 to about 500 fs can be used. Once the pulses are compensated 
(transform-limited), the laser can be focused to produce peak intensities from 
about 10 12 to about 10 18 W/cm 2 , depending on the input energy. 

[0045] This single beam method is capable of retrieving the magnitude 
and sign of second and third order phase modulation (in other words, linear and 
quadratic chirp) directly, without iteration or inversion procedures. MIIPS 
achieves accurate phase retrieval from chirped ultrashort pulses. For MIIPS, no 
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synchronous autocorrelation, beam splitting, or time delays are required because 
the second harmonic spectrum depends on the relative phases of all frequencies 
within the pulse. The amplitude of the pulse is obtained directly from a 
spectrometer in a communications receiver. In order to precisely determine of 
the phase of all frequency components in a pulse from a fs laser 123 (see Figure 
1), a pulse shaper, such as the one described in A. M. Weiner, "Femtosecond 
pulse shaping using spatial light modulators," Rev. Sci. Instrum. 71, pp. 1929- 
1960 (2000), is employed to introduce a reference phase function designed to 
yield this information directly, as further described hereinafter. The shaped 
pulses are frequency doubled by a thin SHG crystal 507 (see Figure 1) and the 
output is directed to spectrometer 503. 

[0046] In addition to laboratory testing and specimen optic distortion 
analysis, the MIIPS system and method employing this single shot construction 
can also be applied to some communication situations in order to add 
considerably more encoded information into each pulse phase to supply 
additional encoding variables. 

[0047] The MIIPS method is based on the principle that second 
harmonic generation, as well as other nonlinear optical processes, depend on the 
phase function </> (w) across the spectrum of the laser pulse. The phase function 
can be expanded in a Taylor series around carrier frequency Q = co - o)o as 
follows: 

M = M +<P'(ab)a+ V2<t>'\(*)G z + / 6 </>'"{ojq)Q 3 +..., 
[1] 
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where the first two terms provide only the relative (common) phase and a time 
delay, respectively. Only the third and higher terms are responsible for phase 
distortion. These higher terms are retrieved in Ml I PS by superimposing a 
reference phase function on the pulse to obtain, 

(f{Q) = acos(yQ-S) + 0(62) 
[2] 

where the first term is the reference phase function introduced by the shaper with 
maximum phase amplitude a, period /and the absolute position in the spectral 
window S. (/{£2) is given by Equation 1 . 

[0048] The maximum SHG signal as a function of Q is obtained when 

d (p{Q)ldG l =0 A p arame t er j n the reference phase function can be varied to 
obtain a plot from which the phase distortions {<p", can be obtained in the 
laser pulse. The maximum signal in a (wavelength, 8) MIIPS trace describes a 
series of lines given by 

<?max = So+ (/Imax -7TC/ab)&b 2 /(7lC) [y- (j)'" l{OC? So)}, 

[3] 

where S ma x is the position where maximum SHG signal is obtained, Sq = arc 
cos[<z>'7(a^)], and /Uax is the position of the maximum SHG signal. 

[0049] A complete data set, from which a phase function can be 
retrieved, consists of a series of spectra obtained as a function of the parameter 
S. The resulting experimental MIIPS trace shown in Figure 16, contains the 
required information to extract Q'W and higher order terms as follows. First the 
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data is fit to a series of lines which follow /lmax(<?max) as expected from Equation 
3. The quadratic phase modulation (responsible for linear chirp) is determined 
directly from the distances x^ and x 2 between the SHG maxima (see Figure 16), 
according to 

0" =ay 2 arc sin[(x r x 2 )/4]. 

[4] 

Note that the magnitude and sign of </>" are obtained directly from the MIIPS 
trace. Furthermore, the accuracy of the measurement can be improved for small 
phase distortion by decreasing the reference phase function parameters ay 2 . 

[0050] The cubic phase modulation (quadratic chirp) is determined by 
the slope AS/Ay that the maximum SHG features make in the X 8 plane. 
Analytically, cubic phase modulation is given by 

</>":= 0.5 ay 2 ml ax? cos[(x r x 2 )/4] {(AS/Ay), - (AS/Ay) 2 }, 
[5] 

where the slopes are measured in nm' 1 (see Figure 16). Higher order phase 
distortions, such as self-phase modulation and quadratic phase components can 
be obtained from the curvature of the line defined by the maximum SHG 
response. The MIIPS can be programmed to find the phase distortions on the 
laser pulses directly by integration and to introduce a compensation phase 
function that eliminates the distortions. This mode of operation can be used to 
find arbitrary phase deformations and yield transform limited pulses, which in a 
MIIPS scan look like straight parallel lines separated by n. The fit to the 
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experimental data shown in Figure 16 - 17C is given by Equation 3, and the 
phase parameters are extracted with Equations 4 and 5. 

[0051] The version of MIIPS illustrated in Figure 15 uses a thin SHG 
crystal 507, spectrometer 503, pulse shaper 129 and a femtosecond laser 123. 
A fs laser pulse is preferred but, for test data disclosed herein, 50 fs pulses from 
a regeneratively amplified Ti:Sapphire laser are employed wherein the pulse 
energy is attenuated down to ~ 5 uJ. For the test data herein, A 0.3 mm pBBO 
type I crystal is used for SHG 507 and the output is attenuated and directed to 
spectrometer 503 with a cooled CCD detector 511. System 121 further has a 
redirecting mirror 513, two quartz cylindrical lenses 515 (200 mm focal length, 
the upstream one for focusing and the downstream one for collimating). For the 
tests, a spatial light modulator was used for pulse shaper 129 consisting of two 
1 28 LCD elements (which can be obtained from CRI Inc. as model number SLM- 
256). For the test, the pulse shaper is carefully calibrated to provide accurate 
phase delays (better than one degree) with no changes to polarization or 
amplitude. The phase distortions used to obtain the data are generated at the 
pulse compressor after regenerative amplification. Referring now to Figures 13 
and 14, self-ultrafast switching is based on pulse phase modulation in pulse 
shaper 505, a thin SHG crystal 507 causing multiphoton intrapulse interference, 
dispersive optics 523, and CCD camera detector 511. The simplicity and 
accuracy of this method make it practical for the evaluation laser pulses close to 
transform limit and for the evaluation of phase distortion from optical elements. 
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[0052] Additional information about Mil and MIIPS is disclosed in U.S. 
Patent Publication No. 2003-0099264 entitled "Laser system using ultrashort 
pulses" invented by Dantus, et al., and PCT Publication No. WO 02061799 
entitled "Laser system using ultrashort pulses" invented by Dantus, et al.; I. 
Pastirk, J.M. Dela Cruz, K.A. Walowicz, V.V. Lozovoy, M. Dantus, "Selective two- 
photon microscopy with shaped femtosecond pulses," Opt. Express 11, 1695- 
1701 (2003); J.M. Dela Cruz, I. Pastirk, V.V. Lozovoy, K.A. Walowicz, M. Dantus, 
"Multiphoton intrapulse interference 3: Probing microscopic chemical 
environments," J. Phys. Chem. A 2004; Lozovoy V.V., Pastirk I., Walowicz K.A., 
Dantus M., "Multiphoton intrapulse interference II: Control of two- and three- 
photon laser induced fluorescence with shaped pulses," J. Chem. Phys. 118 (7): 
3187-3196 (Feb. 15 2003); and Walowicz K.A., Pastirk I., Lozovoy V.V., Dantus 
M.," Multiphoton intrapulse interference 1: Control of multiphoton processes in 
condensed phases," J. Phys. Chem. A 106 (41): 9369-9373 (Oct. 17 2002); all of 
the patent applications of which are incorporated by reference herein. 

[0053] Multiphoton intrapulse interference is not just about focusing the 
energy. The goal is to determine the field that the molecules experience. The 
control of nonlinear optical processes, using multiphoton intrapulse interference 
can be applied in diverse fields such as photochemistry, communications, and 
medicine. 

[0054] The present invention provides a system and method to 
characterize the spectral phase of femtosecond pulses. This single beam 
method is capable of retrieving the magnitude and sign of linear and quadratic 
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chirp with high resolution. Pulse retrieval is based on analytical expressions that 
yield the phase distortion, without iteration or inversion procedures. Linear and 
quadratic chirp values, and to some extent cubic chirp values, are important 
because there are knobs on the laser that can be used to correct for this 
distortion by mechanically adjusting the grating spacing in the laser beam 
amplifier compressor. The method can be used with very short pulses. This 
adjustment can be automatically controlled with the computer controlled software 
as disclosed in Figure 6a. The method is very versatile, and can be used with 
high or very low intensity pulses for any wavelength for which low cost, off-the- 
shelf SHG crystals exist. MIIPS can also be used by obtaining third or higher 
order harmonics in gases. The maximum signal will also agree with Equation 35, 
making the method useful for the characterization of pulses in wavelength 
regions for which SHG crystals are not available. In summary, uses of Mil and 
MIIPS are as follows: 

Mil can be used to make self-switching pulses as long as they undergo 

one non-linear optical process, such as SHG, sum frequency generation, 

difference frequency generation or four-wave mixing; 

MIIPS can be used to allow automated laser optimization, specifically 

quadratic and cubic phase distortions; 

MIIPS can be used for pulse characterization of arbitrary phase 

distortions; 
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MIIPS can be used to measure the phase modulation induced by optical 
elements and similarly it can be used to measure the thickness of a 
substrate; 

MIIPS can be used for decoding information (address and/or message) 
stored in the phase; 

Shapers operating to optimize the MM phenomenon can encode self- 
decoding messages; 

Mil can be used to prevent three photon damage of DNA from fs pulses; 
and 

Mil can be used to optimize the selective activation of PDT agents 

specifically at a particular depth, avoiding collateral damage. 

[0055] Figure 6b is a flow chart showing an automated pulse chirp 
determination for arbitrary smooth phase distortions. This method is based on 
the use of a pulse shaper and obtaining a PhaseScan, wherein the spectrum of 
the SHG as a function of phase parameter Sfor 0(co)=aCos(Y(o+8). This method 
is non-iterative and it directly obtains the desired values without evolutionary 
learning calculations. Therefore this method is very stable. This method does 
not depend on overlap between two pulses in space and time. Moreover, the 
pulse analyzes itself in a thin SHG crystal. 

Second - Harmonic Generation with Powders 

[0056] Chemical powders, adhered onto transparent quartz carriers, 
are employed in place of thin SHG crystals. The powder embodiment is 
presently preferred to significantly reduce cost in high energy (for example, one 
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nanojoule or greater) applications such as for MIIPS, nonlinear optical 
characterizations and FROG analysis. The chemical powder is preferably 
Potassium Dihydrogen Phosphate (KDP or KD*P) or Beta Barium Borate, which 
is glued onto a glass, microscope slide using a silicone-rubber or cyanoacrylate 
(C5H5NO2) adhesive. Type I phase matching is used. The power particle size is 
preferably between about 0.5 to 20 microns, depending on the focusing length, 
intensity of the laser and the sensitivity of the detector. For microscopy, even 
smaller particles can be employed as long as a very sensitive detector is used. 
Furthermore, the powder approach is advantageous by having large variety of 
random crystal orientations therein which creates improved average results for 
laser pulses shorter than 20 fs. 

Pulse Shaping System 
[0057] The first preferred embodiment of a laser system 21 using 
ultrashort laser pulses of the present invention is generally shown in Figure 1 . 
System 21 includes a femtosecond laser 23, an upstream grating 25, an 
upstream concave mirror 27, a laser beam pulse shaper 29, a downstream 
concave mirror 31, a downstream grating 33, a detection device 35, and a 
personal computer 37. Personal computer 37 has a microprocessor based 
electrical control system, memory, an output screen, a data storage device, an 
input keyboard, and a removable disk. More specifically, the detection device is 
a spectrometer 39. Bursts or pulses of a laser beam 43 are emitted from laser 
23, through the optics 25, 27, 31 and 33, as well as through pulse shaper 29 for 
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detection and sensing by spectrometer 39 for further evaluation, analysis, 
comparison and subsequent control by personal computer 37. 

[0058] The laser is preferably an ultra-short femtosecond laser that can 
deliver high peak intensity (with a typical peak greater than 10 10 watts/cm 2 ) which 
preferably emits laser beam pulses of less than 100 femtosecond duration, and 
more preferably at or less than 50 femtoseconds, and for certain applications 
even more preferably at or less than 10 femtosecond duration, for each pulse 
burst or shot. The intense optical pulses that are required to modify material are 
formed in a Kerr-Lens mode locked titanium sapphire oscillator. Such lasers are 
capable of producing hundreds of nanometers of coherent bandwidth, although 
only about 50 nm are typically used. The output may be amplified in a 1 kHz 
regenerative chirped pulsed amplifier. The output pulse is typically 100 fs long 
with a central wavelength of 800 nm and total pulse energy of 0.1 to 1 mJ. 
Preferred lasers include: the Kapteyn and Murnane femtosecond laser oscillator, 
which can produce less than 15 fs pulses at 100 MHz; the Hurricane model from 
Spectra Physics Inc., which is diode pumped and gives 0.8 mJ per pulse with 
sub-50 fs pulses at 1 kHz; and the CPA-2001+ model from Clark-MXR Inc., 
which gives 1.3 mJ per pulse with sub-150 fs pulses at 1 kHz, pumping a Clark- 
MXR Inc. non-collinear parametric amplifier (hereinafter "NOPA") which produces 
0.2 mJ per pulse, and is capable of generating sub-20 fs pulses. This NOPA 
system can even produce pulses between 10 fs and 4.5 fs. 

[0059] A Fourier plane pulse shaper is preferably used with the present 
invention for the transmissive construction illustrated with this embodiment. 
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Ultra-fast laser pulses contain from one to fifty optical cycles, and last only a few 
femtoseconds. This is much faster than most current electronics and therefore 
shaping with fast time gates is very difficult. On the other hand, because of the 
uncertainty principle, the optical spectrum spans tens to hundreds of 
nanometers. Such a large bandwidth is relatively easy to measure and to filter, 
and there are several techniques to shape the spectrum in the frequency domain, 
and thereby shape the temporal pulse upon recompression. 

[0060] In order to access the frequency domain and the individual 
frequency components that comprise the pulse, a geometric arrangement is 
employed, using two back-to-back spectrometers. The spectrometers are 
especially designed to introduce no net temporal dispersion: that is, all colors 
pass through the spectrometers within the same amount of time. The first 
spectrometer (including grating 25 and mirror 27) spreads the unshaped pulse 
spectrum along a line according to its dispersion function y(a). The light 
intercepts spatial amplitude and phase mask pulse shaper 29 at this point. The 
mask output then forms the entrance to a second spectrometer (including grating 
33 and mirror 31) which recombines the colors into a single shaped pulse. 

[0061] The heart of pulse shaper 29 is the programmable 256 pixel 
liquid-crystal mask (consisting of two overlapping 128 pixel liquid crystal arrays) 
that is placed at the Fourier plane. For the applications envisioned herein, the 
mask must be capable of shifting the phase of individual frequencies. For 
alternate embodiment pulse shapers, a different electronically programmable 
mask that is capable of controlling phase has been demonstrated: a liquid 
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crystal display (hereinafter "LCD"), an acousto-optic modulator (hereinafter 
"AOM"), a deformable mirror, and a permanently deformed mirror. A LCD pulse 
shaper can be obtained from CRI Co. and has a modulator electronic driver. 

[0062] The AOM consists of an anti-reflection coated Tellurium Dioxide 
(Te0 2 ) crystal with a piezo electric transducer glued onto one end. The central 
frequency of the acoustic wave is ac/2rc = 200 MHz. The acoustic velocity vs in 
the crystal is 4.2 km/s and the light pulse spends less than 10 ps in the crystal, 
so the acoustic wave moves less than 0.002 X acoustic during the transit of the 
light field through the crystal. Since the acoustic wave is essentially frozen as 
the optical pulse travels through the crystal, the complex amplitude of the 
acoustic wave traveling through the crystal in the y direction, A(t)cosact = 
A(y/vs)cosact, is mapped onto the optical field E(a) as it passes through the 
AOM. If some of the dispersed optical field encounters a weak acoustic wave, 
that frequency is attenuated; if the acoustic wave carrier is shifted by phase 
angle 0, that phase shift is imposed on the optical field. This pulse shaper has a 
total efficiency of about 20% including the diffraction efficiency of the AOM and 
the diffraction efficiency of the gratings. The diffracted light is used and the 
undiffracted "zero order" beam is blocked, to allow full modulation of both 
amplitude and phase in the shaped beam. The shaped beam than has the form 
E shaped (to) = E inpul mxa&W™' [6] 

where a((o)e mw) = A[y((o)/v s ];a is the frequency, and e is a 
constant. 
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Fixed pulse shaping optics, such as chirped mirrors, can also be employed as 
will be discussed further hereinafter with regard to Figure 23. 

[0063] The transform-limited pulses (hereinafter TL"), having all their 
frequencies in phase, are fed into the pulse shaper where curved mirror 27 
focuses the spectrum onto Fourier plane 29. Changes in the phase o and 
amplitude A of the spectral components indicated by the computer are used to 
tailor the laser pulse before reconstruction with second curved mirror 31 and 
grating 33. Once compressed, the shaped pulse is directed to spectrometer 39 
for evaluation. The Fourier transform relationship between the time and the 
frequency domain allows us to calculate the necessary mask to create a certain 
shaped pulse. These calculations are based on 

fM = T~l fit)e' lma dt 

and 

/«=£ f(v)e- i2ma dv 

where v is the frequency in wave numbers, t is the time, and c is the speed of 
light. 

[0064] In this embodiment, the phase and amplitude masks of the 
pulse shaper are controlled by the computer wherein the laser pulse shape takes 
a dynamic role. The microprocessor within personal computer 37 will then 
control laser 23, receive an essentially real time feedback input signal from 
spectrometer 39, and then perform calculations, comparisons and evaluations, 
and possibly automatic variation of subsequent pulse shapes. These automated 
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steps can be substituted with manual user calculations and decisions if desired 
based on personal computer outputs. 

[0065] As applied to all of the applications herein, selective control of 
one and multiphoton processes in large molecules, including proteins, is possible 
using a simple pulse shaping method that is based on taking maximum 
advantage of the multiphoton intrapulse interference caused in short pulses with 
large bandwidths. The results show an extraordinary level of control that is 
robust and sample independent, with contrast ratios near two orders of 
magnitude (clearly visible with the naked eye). Such large contrast ratios allow 
for more precise cancellation control of undesired photons and other laser beam 
characteristics, such that nonlinear transitions induced by each pulse are 
controlled. Because simple phase functions can be incorporated into a passive 
optical component such as mirror 129 (see Figure 5), these applications do not 
require the complexity and expense of computer controlled pulse shapers after 
initial set up, although systems can still be employed. 

[0066] A fs-pulse shaper arrangement can be used to achieve 
background free functional imaging (pH, Na or Ca concentration gradients, 
electric fields, charge, fluorescent probes, nanoclusters, or quantum dots, 
chemical composition) by taking advantage of the selective excitation afforded by 
multiphoton intrapulse interference. For example, A. H. Buist, et al., "Probing 
microscopic chemical environments with high-intensity chirped pulses," Optics 
Letters 24, 244-246 (1999). Buist, et al., showed that linear chirp can be used to 
crudely distinguish the pH environment of a pH-sensitive dye. With phase 
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modulation, and specifically taking advantage of multiphoton intrapulse 
interference, the present invention can achieve much more sensitive pH 
sensitivity with a greater number of pH-sensitive dyes. Using the same principle, 
dyes that are sensitive to sodium, calcium or other chemical gradients including 
also charge and can also be probed selectively. Alternatively, multiple probes 
such as dyes, nanoclusters or quantum dots can be selectively excited through 
two- or three-photon excitation. 

Binary Phase Shaping 

[0067] In the low intensity regime, laser control is dominated by 
interference of different nonlinear optical pathways connecting the initial and final 
states. The challenge is finding the proper phase for each frequency within the 
pulse to achieve constructive interference at the desired pathway and destructive 
interference elsewhere. Consider two-photon excitation of fluorescent probes, as 
used in two-photon microscopy, as the target for optimization and envision two 
chromophores with different two-photon absorption spectra. The goal is to 
achieve selective excitation by "focusing" the energy available at a specific region 
of the two-photon spectrum, while minimizing the energy outside of the desired 
spectral window as shown in Figure 8. 

[0068] The broad bandwidth second harmonic spectrum from 
transform-limited pulses is represented by a Gaussian curve (thin line). The 
objective is to introduce phase modulation to cause the two-photon spectrum to 
be intense only inside the window defined by frequency 2&b and width W, and to 
minimize the background B outside the window. The contrast ratio C is defined 
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as the intensity inside the window divided by the intensity of light outside the 
window 

[0069] The phase between photons of different frequencies takes only 
two values preferably 0 or preferably n to maximize or minimize a given pathway. 
Any two values whose difference is n work equivalent^ well. Based on the 
theoretical and experimental confirmation of this, the method is defined as binary 
phase shaping (BPS). In a preferred embodiment, BPS is used to solve the 
problem of selective multiphoton excitation with ultrashort laser pulses. 

[0070] In order to control two-photon excitation, one needs to control 
the non-linear power spectrum of the laser & 2) {co), which can be measured by 
obtaining the second harmonic spectrum generated using a thin second 
harmonic generation (SHG) crystal. Selective two-photon excitation is possible 
when one is able to tune the narrowed non-linear power spectrum to optimize 
excitation of one chromophore versus another. 

[0071] The effect of spectral phase modulation on SHG has been 
studied by a number of groups and can be divided into broad and narrow phase 
matching bandwidth. Sinusoidal phase modulation, a common function used for 
laser control, cannot produce contrast ratios greater than 0.5, and as the window 
is tuned away from the central frequency the contrast drops below 0.1 . 

[0072] Consider two phases, 0 and n, then the symmetry becomes 
clear. To maximize the SHG intensity at a frequency 2co c ,. the spectral phase 
needs to be symmetric or antisymmetric about coc, such that the frequencies 
interfere constructively. To minimize the background intensity at all other 
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frequencies, the spectral phase must be asymmetric with respect to all other 
frequencies away from co c , so that destructive interference is maximized. As 
described by others, binary phase filters masks or deformable mirrors or SLMs 
produce waveforms with symmetrical intensities only. (See A.M. Wiener, 
"Femtosecond pulse shaping using spatial light modulators," Rev. Sci. Instrum., 
Vol. 71(5), page 1934.) In another preferred embodiment, prime numbers are 
used to generate the quasi-random phase changes required. The mask that is 
used to modulate the pulses is shown in Figure 9, and is designed for a 128-pixel 
modulator. As shown in Figure 9, a phase mask proposed based on the 
symmetry requirements of the problem, using the quasi-randomness of prime 
numbers. This mask is reflected about pixel 64 to obtain the values of pixels 65- 
128, and is designed to obtain a narrow second harmonic signal at the center of 
the spectrum, for example. Other preferred embodiments employ a 128-pixel 
SLM and still other embodiments employ a SLM with greater than 256 pixels. 

[0073] A titanium-sapphire oscillator which can be obtained from K&M 
Labs laser system capable of generating pulses as short as 10 fs after a double 
pass prism compressor is preferably employed. The spectral phase of the pulse 
is tailored using a computer-controlled pulse shaper. Preferably, the pulses are 
centered near 800 nm. The spectral phase of each pulse was corrected using the 
Mil phase-scan (MIIPS) method, which compensates phase distortions to obtain 
transform-limited (TL) pulses. The binary phase was introduced as an addition to 
the compensation phase. 
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[0074] The shaped laser pulses, with energy ~ 0.5 nJ per pulse and 87 
MHz repetition rate, are focused mildly, to a spot size of ~ 100 microns in 
diameter, on a 20 micron thin beta barium borate ((3BBO) type I SHG crystal. The 
frequency-doubled light was collected with an optical fiber and dispersed on a 
compact spectrometer, preferably obtainable from Ocean Optics. 

[0075] Before introducing phase modulation, spectral amplitude 
restriction was considered and the spectrum of the laser was narrowed using a 
slit at the Fourier plane. The resulting spectra are shown in Figure 10a. In Figure 
10b the SHG intensity obtained after amplitude restriction is shown and note that 
the energy inside the window depends on the spectral width squared. For a 
spectral width of 1 0% of the available bandwidth, amplitude restriction produces 
100 times less second harmonic intensity than TL pulses. Although the contrast 
ratio for this case is favorable, the low intensity obtained at the desired 
wavelength makes this method experimentally impractical. 

[0076] The phase mask shown in Figure 9 is programmed on the SLM 
and it shows dramatic narrowing of the SHG spectrum as illustrated in Figure 4. 
The contrast ratio as defined in Figure 1 for this mask is 2.5 when the SHG peak 
is centered. By shifting the position of the mask on the SLM, hence tuning the 
center of symmetry, and thus tuning the SHG spectrum as shown in Figure 1 1 . 

[0077] The dashed line in Figure 11 corresponds to theoretical 
calculations. The absolute value of the spectral amplitude \E{cd)\ of the electric 
field is calculated from the experimental power spectrum of the fundamental 
pulse l{co) using \E(co)\ = /(a>) 05 . To simulate the experimental results a double 
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Fourier Transform method is used. The electric field in the time domain E(f) is 
calculated as the Fourier image of the complex spectral amplitude in the spectral 
domain, with the formula 

E(t) = J| E(co)\exp[i(p{Q))]exp(-icot)da), [8] 
where the spectral phase <fl(co) is the function that is introduced by the SLM. The 
power spectrum of the SHG is calculated using 

/shgM = |jE(f) 2 exp(/Vyf)df . [9] 
The SHG amplitude is normalized using the maximum of the SHG intensity 
calculated for TL pulses, (j)(co) = 0. 

[0078] A simple evolutionary learning computer program (ELC) is 
programmed and it assumed a Gaussian electric field \E{co)\ corresponding to a 
10 fs pulse centered at 800 nm. Eight points are used to represent each pixel in 
order to simulate more closely the experimental resolution of our setup. The 
second harmonic intensity is calculated according to 

Ishg(2o) c ) = \lE(Q) c -Q))E{o) c +Q))6cif [10] 
normalized to the maximum SHG amplitude for TL pulses. 

[0079] BPS simplifies the calculations, especially if it is assumed the 
amplitude of the electric field to be a constant, that is, the spectral power is set 
equal to 1 in the allowed spectral region. Each spectral component of the electric 
field, linearly dispersed in the frequency domain, can be represented as a binary 
value (± 1) determined by b k = exp(i<t> k ), for <h< = 0 or n respectively. The intensity 
of the SHG signal measured at frequency 2w k can be calculated with the formula 
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7=0 



[11] 



where the integral in Equation 10 is now replaced by a discrete sum, N is a 
parameter that depends on details of the model such as number of pixels. The 
problem of spectral selectivity can now be formulated as finding a vector l\ such 
that S k = 1 for cok = cue and S k is minimized at all other frequencies. 

[0080] In Figure 12, theoretical results from the phase described in 
Figure 9 (black) are compared to the results of the ELC (dashed) which was 
initiated with the proposed phase. The evolutionary calculations used 100 
individuals and one parent through 30 generations using single double and triple 
symmetric bit flips. The result from the ELC gives -2.5 greater contrast. Using an 
SLM with more pixels, and grouping those pixels such that each group covers the 
same frequency range, improves the implementation and result in greater 
contrast, an observation based on our calculations. The goal of narrowing the 
SHG spectrum was motivated by the possibility of selective two-photon 
microscopy. To achieve this goal two conditions had to be satisfied, maximizing 
the nonlinear power spectrum at 2&b, and minimizing it elsewhere. Maximization 
has been investigated previously by Silberberg and coworkers, who identified the 
condition § (ftb-<y) = -<|> (ah+oj) for continuous functions. (See: N. Dudovich, B. 
Dayan, S.M.G. Faeder, Y. Silberberg, "Transform-limited pulses are not optimal 
for resonant multiphoton transitions," Phys. Rev. Lett. 86, 47-50 (2001) and D. 
Meshulach, Y. Silberberg, "Coherent quantum control of multiphoton transitions 
by shaped ultrashort optical pulses," Phys. Rev. A 60, 1287-1292 (1999)). Zheng 
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and Weiner explored the narrow phase matching bandwidth SHG output at a 
single frequency, a problem that is analogous to two-photon excitation of an 
atom, using binary encoded pulses as used in communications. They found that 
limiting the phase to values of zero and n led to the condition 0 (cuk-co) = <J> {ok+a)). 
(See: Z. Zheng, A.M. Weiner, "Spectral phase correlation of coded femtosecond 
pulses by second-harmonic generation in thick nonlinear crystals," Opt. Lett. 25, 
984-986 (2000) and Z. Zheng, A.M. Weiner, "Coherent control of second 
harmonic generation using spectrally phase coded femtosecond waveforms," 
Chem. Phys. 267, 161-171 (2001)). Under both shaping conditions the 
maximum signal approached that is obtained with TL pulses. 

[0081] The problem of spectral narrowing by pulse shaping under 
broad phase matching SHG conditions, or selective excitation of large organic 
dyes and chromophores requires minimization of the nonlinear power spectrum 
away from 2(o c . As expressed in Equation 11, there is an analogy with 
convolution. In this case, some binary functions have the property of giving a 
sharp 'delta function' convolution signal in one location and very low background 
elsewhere. 

[0082] In principle, the solutions found are members of the set of 
solutions that could be obtained by arbitrary phase and amplitude pulse shaping. 
For a pulse shaper with N pixels, one can generate (P *A) N shaped pulses, 
where P and A are the number of different phases and amplitudes a pixel can 
take. If it is assumed 100 pixels, each taking 10 different amplitude values and 
100 different phase values, the number of different pulses is of order of 
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magnitude 10 300 . This number is extremely large, therefore, while in principle, the 
field exists to achieve the desired photonic transformation or excitation, finding it 
is a great challenge. 

[0083] The periodic nature of electromagnetic waves results in a great 
deal of redundancy in pulse shaping because nonlinear optical processes do not 
depend on the absolute phase or a linear variation of the spectral phase. This 
equivalence is expressed by §ctK=>§{a))+a+bco, where a and b are constants. This 
redundancy is filtered out by programming an ELC that works on the second 
derivative of the phase. The actual phase that is used in the SLM is obtained by 
integration setting a = b = 0. An ELC was used to optimize smooth phase 
functions for spectral narrowing, but could not obtain a contrast ratio greater than 
unity. 

[0084] The advantage of BPS is that computational redundancies are 
greatly reduced. For BPS and 128 active pixels the search space is reduced to 
2 128 . If there is two-fold symmetry, for example two-photon excitation, then the 
number is reduced to 2 64 . The final search space is of size 10 19 , a number that is 
at least 281 orders of magnitude smaller than would be considered for arbitrary 
phase and amplitude pulse shaping as discussed above. The resulting space is 
small enough that all possible outcomes could be computed, and a large portion 
evaluated experimentally. A simple ELC, such as the one used here can quickly 
converge towards significantly improved solutions. 

[0085] BPS may have significant technological advantages. A 
retardation equivalent to n is easy and fast to obtain and calibrate. Permanently 
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etched masks can be made in advance and used for specific applications such 
as selective two-photon microscopy. Scanning the mask can yield two-photon 
excitation spectra. 

[0086] It has been demonstrated that laser control, especially with two- 
photon transitions, can be addressed with binary phase shaping. BPS makes it 
simple to analyze the problem and to propose rational solutions, as 
demonstrated here with a phase mask fabricated by the quasi-random gaps 
between prime numbers. Thus, a simple ELC is used to improve on the 
proposed solution efficiently because of the greatly reduced search space. 
Applications Using Binary Phase Shaping 

[0087] Nonlinear optical processes, having a high-order dependence 
on light intensity, play a central role in a number of key technologies, among 
them multiphoton microscopy, multiphoton photodynamic therapy, multiphoton 
microlithography, optical switching, chemical sensors, and selective 
photochemistry. A sizable research effort is underway to improve the efficiency of 
these processes, focusing primarily on materials with large nonlinear optical 
susceptibility. Being able to enhance the efficiency for nonlinear excitation, or, in 
some cases, being able to suppress nonlinear optical distortion, remains a top 
priority in the design of photonic technologies. Research groups have been 
using laser systems with characterization and shaping technology to control 
nonlinear optical processes in large molecules, including proteins. Using the 
amplitude of the pulse and the spectral phase across its bandwidth, the power 
spectrum of the field can be determined. The spectral phase introduced by the 
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shaper controls the way frequencies in the pulse add up to produce the higher 
nonlinear frequency components. The phase-dependent probability is used to 
manage the frequencies over which nonlinear optical processes are induced by a 
laser pulse. 

[0088] Phase modulation introduces changes in the power spectrum of 
multiphoton processes. If the phase function <D in introduced in the SLM and 
scanned 6 from 0 to 2n, the power spectrum responsible for two-photon induced 
fluorescence changes accordingly. This approach is usable to achieve selective 
multiphoton microscopy. The system produces selective excitation of 
microscopic polystyrene beads doped with blue or green fluorescent dyes. 
Transform-limited pulses excite multiphoton transitions from both types of beads; 
however, the spectral phase of the laser pulse is modulated such that the pulse 
excited only specific types of molecules, providing selectivity. For these 
experiments, the output from the pulse shaper is focused directly on the sample 
being imaged by a microscope; the spectrum and intensity of the incident pulses 
remains unchanged in all cases. 

[0089] Applications for this method include all the techniques that use 
multiphoton excitation. This method is used for pulse characterization and 
compensation, selective control of two and three-photon excitation, selective 
microscopic environment sensing and imaging, and selective multiphoton 
microscopy and imaging. MIIPS and BPS offer a powerful method for pulse 
characterization and compensation of ultrashort laser pulses. All applications 
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involving nonlinear optical processes will benefit from the control that spectral 
phase manipulation brings through multiphoton intrapulse interference. 

Multi-Photon Microscopy 
[0090] Two-photon microscopy provides significant possibilities for 
fluorescence imaging and photochemistry. It offers attractive advantages, 
including higher resolution, background-free signal, lower background scattering, 
better penetration in thick samples, and reduced photon-induced damage, which 
arise from the basic physical principle that the absorption depends on the square 
of the excitation intensity. Two-photon microscopy is amenable to multiple-probe 
staining, whereby two-photon transitions excite different probe molecules that 
emit at different wavelengths, and for functional imaging of living cells. See U.S. 
Patent No. 5,034,613 which issued to Denk, on July 23, 1991; U.S. Patent No. 
6,166,385 which issued to Webb, on December 26, 2000; U.S. Patent No. 
6,344,653 which issued to Webb, on February 5, 2002; U.S. Patent No. 
5,759,767 which issued to Lakowicz, on June 2, 1998; and W.R. Zipfel, et al., 
"Nonlinear magic: multiphoton microscopy in the biosciences," Nature 
Biotechnology, 121 (11): 1369-1377 (Nov. 2003); all the above patents are 
incorporated herein by reference. Phase-modulated femtosecond pulses can 
selectively excite one type of probe molecule only, leaving the others in their 
ground state. Multiphoton excitation is achieved by multiphoton intrapulse 
interference (Mil) and this can be accomplished efficiently using binary phase 
shaping.. Selective excitation is used to enhance contrast and achieves 
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functional imaging of samples stained with fluorescent probes sensitive to their 
microscopic chemical environment. 

[0091] Preferably, a titanium-sapphire oscillator, which can be obtained 
from K&M Labs, is employed in the present invention laser system and is 
capable of generating pulses as short as 10 fs after a double pass prism 
compressor. The amplitude profile and spectral phase of the pulse is tailored 
using a computer-controlled pulse shaper of a design as described above. The 
shaper uses a pair of SF-10 prisms as the dispersion elements and 0.4 m 
cylindrical mirrors for collimation. A schematic of the experiment is shown in 
Figure 7. Phase modulation is introduced by a computer-controlled spatial light 
modulator (SLM-256, CRI Inc.), located at the Fourier plane of the shaper setup. 
The pulses are centered near 800 nm. The spectral phase of each pulse is 
calibrated using the multiphoton intrapulse interference phase-scan method 
(MIIPS), first to obtain transform-limited (TL) pulses and then to introduce the 
desired phase. The spectral phase function used for this experiment is given by 

(p(Q)=acos(yQ-5) [12] 
where Q is the frequency detuning measured in s" 1 from the center wavelength of 
the laser pulse. The shaped laser pulses, with an energy of 0.2-1 nj per pulse 
and 87 MHz repetition rate, are focused mildly, to a spot size of - 200 microns in 
diameter, on the different samples. Two-photon induced fluorescence is 
collected by an infinity corrected apochromatic microscope objective (APO 100X, 
Mitutoyo) and imaged into a liquid nitrogen cooled CCD detector using an 
InfiniTube, preferably a Proximity Series model obtained from Infinity Photo- 
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Optical. Experiments were carried out on three different types of samples. The 
first sample is prepared with acrylonitrile-vinylidene chloride polymer cured in 
acidic (pH 6) or basic (pH 10) buffered conditions and stained with the pH- 
sensitive fluorescent probe 8-hydroxypyrene-1,3,6-trisulfonic acid (HPTS) (10" 3 M) 
(Sigma-Aldrich). The second sample consists of pieces of 
polymethylmethacrylate (PMMA) doped with C540 and R6G fluorescent probes 
(10' 4 M). The third sample consists of a mixture of 10 urn blue and 15 urn green 
fluorescent FluoSpheres (Molecular Probes). The images are shown in black and 
white without optical filters in the collection optics or post image manipulation. 
The background, detector counts when the laser is off, is subtracted. 

[0092] Figure 19 shows the experimentally measured and theoretically 
predicted dependence of the contrast ratio of intensity of two-photon laser 
induced fluorescence between samples with pH 1 0 and pH 6 for the solutions of 
pH sensitive dye in order to demonstrate the Mil principle. Figure 19 shows the 
experimental results and theoretical predictions of selective two-photon excitation 
of HTPS solutions at different pH: (a) the laser spectrum plus a phase function; 
(b) calculated power spectrum of E 2 (t) for TL pulse; absorption spectra of HTPS 
at pH 10 (right side) and pH 6 (left side). 

[0093] 10' 5 M aqueous solutions of HPTS salt in different pH buffers 
are the exemplary samples. The laser spectrum remains unchanged throughout 
the experiments as spectral phase is modulated. The excitation spectra of HTPS 
exhibits large pH dependent differences while the emission spectrum centered at 
515 nm is independent of pH. These measurements are obtained using 23 fs 



38 



6550-000057/CPE 

transform-limited pulses that were shaped to have a Gaussian spectrum 
centered at 842.6 nm with a bandwidth of 48 nm and spectral phase defined by 
Eq. [12] with a = 1.5tt, y = 20 fs. Spectral intensity |E(Q) 12 and phase <t>(Q) of 
laser pulse are presented in Figure 19a. Figure 19b shows the calculated 
spectrum of the electric field squared E (2) (A) = E(A/2+Q)E(A/2-Q)dQ [13] where 
Q=2kc(1/A-1/A 0 ) is the spectral detuning from the carrier frequency 2tcc/A0 of the 
laser pulse and E(Q) = |E(Q) |ei4,(n) . Simulations of the experimental signal are 
performed using the following equation with no adjustable parameters: 

S< 2 >aJg (2) |E< 2 >(A)| 2 dA [14] 
where y (2) (A) is two-photon absorption spectra of HTPS at different pH (single 
photon spectra are shown in Figure 19b. The spectrum of the pulse |E(Q)| 2 and 
the absorption spectra g(A) of HPTS are measured experimentally. 
Intramolecular dynamics, such as molecular vibrations, are neglected because 
they have a minimum contribution to Mil. Because HPTS does not have 
inversion symmetry, its one-photon absorption spectrum g(A) is expected to be 
similar to its two-photon absorption spectrum y (2) (A). The contrast ratio between 
the measured intensity of fluorescence at pH 10 and pH 6 is presented as a 
function of 6. The experimental data obtained by excitation of the dye solutions 
in methanol (10 3 M in water) generally agree with the calculated response. 

[0094] Mil modulates the amplitude of the power spectrum of E 2 (t) as 
the spectral phase of the laser pulse <t>(A) is scanned. Note that the maximum 
amplitude in the power spectrum in Figure 19d shifts from longer to shorter 
wavelengths, giving the possibility for selectively exciting the fluorescent probe in 
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different pH environments. A solid sample of HTPS is prepared with different pH 
regions and imaged using a microscope. Figure 20a shows an image obtained 
with TL pulses. The spectrum of the pulse and the spectral phase (flat line) are 
shown in the right panels. Figure 20b shows an image obtained with laser pulses 
optimally shaped for selective excitation of HPTS under acidic environment. The 
spectral phase used to acquire this image is shown in the middle panel. 

[0095] The experimental results in Figure 21 show selective excitation 
of different fluorescent probes. Figure 21a shows the data obtained using TL 
pulses. The two pieces of PMMA being imaged show similar amounts of two- 
photon laser induced fluorescence. When the spectral phase of the pulse is 
modified to optimize selective excitation of the C540-doped sample, only the top 
fragment fluoresces were observed, as shown in Figure 21 d. When the spectral 
phase of the pulse is optimized for selective excitation of the R6G-doped sample, 
Figure 21c, strong fluorescence was observed from the bottom fragment while 
the other side, containing the C540 probe, does not fluoresce. Note that the 
spectrum and intensity of the laser remain constant during the experiment for all 
three cases. The only change is the subtle modulation of the spectral phase of 
the pulse, which controls multiphoton intrapulse interference. Contrast ratio 
(IR6G/IC540) (normalized intensities corrected to background) ranges from 1:9 
for pulses optimized for C540 to 8:1 for pulses optimized for R6G. Selective 
multiphoton excitation minimizes the cross talk between different fluorescent 
probes and does not require the use of multiple filters and detectors for imaging. 
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[0096] The previous experiment shows selective two-photon excitation 
of different probe molecules. To further illustrate the potential of coherent control 
for functional imaging, selective excitation of microspheres is shown which is 
typically used for targeted staining of biological samples. Fluorescence from two 
different microspheres is shown in Figure 22. The top image shows the green- 
fluorescent microsphere which has an absorption maximum at 450 nm, while the 
bottom image represents the blue-fluorescent microsphere, which has an 
absorption maximum at 365 nm. The image in Figure 22a is acquired with TL 
pulses. Under these conditions, large two-photon induced fluorescence signal is 
observed from both spheres. For the image in Figure 22d, the phase of the 
pulses is modulated such that the intensity of fluorescence from the blue 
microsphere is maximized and that of the green one minimized. 

[0097] For the image in Figure 22c, the phase of the pulses is 
modulated such that the fluorescence from the green microsphere is maximized 
and that of the blue one minimized. The observed contrast ratios (Iblue/Iqreen) 
are 1:3 and 4:1 for the different phases. 

[0098] Selective excitation with significant contrast ratios has been 
achieved here by optimizing the overlap between the power spectrum of E 2 (t) 
and the two-photon absorption spectrum. The addition of a computer-controlled 
pulse shaper to the multiphoton microscope provides a number of important 
advantages. First, the pulse shaper is used to compensate unwanted phase 
distortions at the sample. Linear chirp, for example, has been shown to reduce 
signal intensity in two-photon microscopy. With a pulse shaper, linear, quadratic, 
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cubic and higher order chirp is compensated to obtain the most efficient 
excitation. Second, the pulse shaper is used to control the output spectrum of the 
laser pulses by amplitude modulation. Third, the pulse shaper is used for 
selective probe excitation as shown here. Because the spectrum of the laser 
remains constant, phase modulation does not affect one photon processes such 
as absorption, reflection and scattering. Selective excitation minimizes possible 
cross talk between different fluorescent probes in the sample. Finally, the pulse 
shaper is used to prevent three-photon and higher order nonlinear optical 
processes such as continuum generation. Higher order processes usually lead to 
sample degradation, and in the case of living samples to DNA damage. 
Suppression of three-photon transitions of four orders of magnitude has been 
achieved using the Mil and BPS methods and this suppression can be coupled 
with optimization of two-photon signal from living specimens. 

[0100] This method can be used to selectively excite either different 
probe molecules or identical probe molecules in different environments. In 
addition, this method can be used for selective excitation of luminescent probes 
such as quantum dots, metallic nanoparticles, and single molecules. For 
example, see, D.R. Larson, et al. " Water soluble quantum dots for multiphoton 
imaging in vivo," Science 300 1434-6, (May 30, 2003); and D.A. VandenBout, et 
al., "Discrete intensity jumps and intramolecular electronic energy transfer in the 
spectroscopy of single conjugated polymer molecules," Science 277, 1074-1077 
(1997). The same principle can be extended to achieve functional imaging of 
semiconductor microchips by two-photon laser induced conductivity. For 



42 



6550-000057/CPE 

examples see: C. Xu, et al., "Two photon optical beam induced current imaging 
throughout backside of integrated circuits," Appl. Phys. Lett. 71, 2578-2580 
(1997) and D.L. Osborn, et al., "Spectral and intensity dependence of spatially 
resolved two- photon conductivity defects on a GaAsP photodiode," J. Appl. 
Phys. 89, 626-633 (2001). Having a pulse shaper for multiphoton microscopy 
provides the flexibility of selective probe excitation or maximum signal 
enhancement by controlled modulation of the spectral phase of the femtosecond 
pulses. Even when for fluorescent labels with very similar absorption spectra, 
pulse shaping is shown capable of selective excitation. This level of selective 
excitation and enhancement can be adapted to different modes of two-photon 
and three-photon microscopy. 

Multi-Photon Photo Polymerization 

[0101] Microlithography involves the use of UV light to initiate 
polymerization. Two-photon induced polymerization has the advantage that it 
permits 3D polymerization of smaller features (down to 100 nm). Photonic band 
gap materials (PBGM) are 3D constructions with features similar to the 
wavelength of light that exhibit very interesting behavior. Two-photon 
microlithography is one convenient method for the preparation of very 
sophisticated PBGMs. 

[0102] Two-photon absorption can be used in three-dimensional 
lithographic microfabrication (3DLM) through two-photon-induced polymerization 
(TPIP). This is due to the fact that simultaneous two-photon absorption requires a 
very high photon flux, which is only present at the point of the focus. Thus the 
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TPI polymerization is confined to the focal volume. This high spatial resolution 
contributes to the ability of TPIP not only to scan the laser in the x and y direction 
but also to change the focal plane (z) without overwriting existing features. 
Therefore, 3DLM is obtained by a single processing step. 3D polymeric 
structures include a photonic band gap structure, waveguide structures and a 
micro-channel structure. 

[0103] Examples of molecules that can be used in 3DLM have the 
general structure D-n-A-n-D, wherein D is a terminal group, n is a phenylimine 
function and A is an electron-accepting cyano group. Since cyano groups are 
used as the side groups, the increase of electron transfer leads to a significant 
increase in the 6 value relative to electron-donating methoxy groups. 
Additionally, the change of the conjugated backbone through the introduction of 
the polar imino groups will allow the tuning of the linear and nonlinear optical 
responses of quadrupolar derivatives as a result of the influence of the 
chromophoric two-photon absorption properties. Both polar groups (-CHLN-) 
will enhance the extent of the charge transfer from the terminal groups to the 
center, correlated with the two cyano-groups. Upon chromophoric excitation, the 
p-conjugated center acts as an electron acceptor, and intramolecular energy and 
electron transfer will occur more readily from the terminal R2N- groups than via 
intermolecular energy and electron transfer. Radical forms can be easily 
generated in the terminal groups. Moreover, the solubility of the chromophore 
improves with the replacement of -CHLCH- by the polar flexible -CHLN- group. 
Due to the fact that they can be synthesized in one-step with a high yield and 
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initiate polymerization of triacrylate monomer at a high level of sensitivity, makes 
application in two-photon-induced polymerization possible. 

[0104] The following references disclose two-photon photopolymer 
initiators: "New Photopolymers based on Two-Photon Absorbing Chromophores 
and Application to Three-Dimensional Microfabrication and Optical Storage," B. 
H. Cumpston, et al., Mat. Res. Soc. Symp. Proa, Vol. 488, "Electrical, Optical, 
and Magnetic Properties of Organic Solid-State Materials IV," (MRS, Warrendale, 
1998) p. 217; and "Two-Photon Polymerisation Initiators for Three-Dimensional 
Optical Data Storage and Microfabrication," B. H. Cumpston, et al., Nature, in 
press. It is envisioned that multiphoton intrapulse interference can be 
advantageously used to enhance this non-linear photopolymerization. 

[0105] The following references are examples multi photon techniques 
for microlithography and microfabrication: Y.M. Lu, et al. "Highly sensitive two- 
photon chromophores applied to three dimensional lithographic microfabrication: 
design, synthesis and characterization towards two-photon absorption cross 
section" J. Mater Chem. 14 (1): 75-80 (2004); K.D. Belfield, et al, "Two-photon 
photoinitiated polymerization," J. Phys. Org. Chem. 13 (12): 837-849 (December 
2000); B.J. Postnikova, et al, "Towards nanoscale three-dimensional fabrication 
using two-photon initiated polymerization and near-field excitation," 
Microelectron. Eng. 69 (2-4): 459-465 (September 2003); H.B. Sun, et al, 'Two- 
photon laser precision microfabrication and its applications to micro-nano devices 
and systems," J. Lightwave Technol. 21 (3): 624-633 (March 2003); U.S. Patent 
6,316,153 to Goodman, issued November 13, 2001; U.S. Patent 6,259,104 to 
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Baer, issued July 10, 2001; and U.S. Patent 6,288,782 to Worster, issued 
September 11, 2001; all of the above patents are incorporated herein by 
reference. 

[0106] Currently, two-photon microlithography and related techniques 
can only initiate one type of polymer. If two different types were needed, one 
would need to rinse, change the monomer mixture, find the position of the feature 
with nanometer accuracy, and make the new feature. 

[0107] In one embodiment, BPS is used in microlithography and 
related techniques such as micromachining and microfabrication to control the 
polymerization of two different polymers. Being able to alternate two different 
polymers allows greater flexibility in the construction of nanometer features such 
as, but not limited to, microelectromechanical systems (MEMS). Controlling the 
deposition of two different polymers is enabled by the control of what wavelength 
two-photon excitation takes place. 

Quantum Computing 

[0108] One such important application of the concepts of coherent 
control is emerging in the area of quantum computing. The selectivity in 
excitation, which is offered by shaped laser techniques, may also serve as the 
building blocks for the development of the first practicable quantum computer. In 
a prior optical scheme for quantum computing, computer controlled pulse 
shaping where information storage and retrieval through quantum phase is 
required. In case of an eight-state Rydberg atom wavepacket, the prior scheme 
can store information as a quantum phase in one or more flipped state, which 
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could be subsequently retrieved in a single step in agreement with the Graver's ; 
this prior scheme is disclosed in J. Ahn, et al., Science 287 (2000), p. 463; and 
N. Bhattacharya, et al., Phys. Rev. Lett. 88 (2002) 137901-1. 

[01 09] A typical visualization of a traditional quantum computer network 
would have nodes consisting of quantum storage devices, where information can 
be stored for very long times either in ground or in some metastable excited 
states of atoms, molecules or ions. The quantum information can be transferred 
from one node of the network to the other using photons. The nodes would carry 
out the required computations and also serve as a storage or memory unit. The 
storage time is limited by decoherence time. Transferring quantum information 
between the two nodes without allowing for decoherence is very difficult. There 
are already some proposals in quantum communication to transmit and 
exchange quantum information between distant users, which includes distribution 
of quantum secure key information for secure communication. Teleportation 
allows an arbitrary unknown quantum state to be conveyed from one distant part 
to another with perfect fidelity by the establishment of a maximal entangled state 
of two distant quantum bits. However, the bottleneck for communication between 
distant users is the scaling of the error probability with the length of the channel 
connecting the users. The error results from amplitude and phase damping. 

[0110] Advantageously, optical pulse shaping of the present invention 
is an attractive route to quantum computing since shaped pulses can be 
transmitted over optical hardware and the same infrastructure can be used for 
computation and optical information transfer. The shaped pulses are split into a 
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number of different parts which can carry different train of pulses at different 
timescales. This provides leverage to control the various nodes where molecular 
systems interact with shaped pulses to carry out various instructions and perform 
quantum computing activity at each node during the pulse duration. This further 
enables the processing of different quantum computational steps at various 
nodes simultaneously, such that the code is transmitted in parallel for distribution 
of the task over the network. At the end of the computation the results are read 
by sending in a "read pulse" and recombining the results. Essentially, this is 
distributed quantum computing over the network using shaped pulses. Currently, 
106 bits can be transmitted/encoded in a single burst of light with the present day 
optical pulse shaping technology. The repetition rate from the laser source is 
about 50 to 100 MHz. Thus, one would be able to use terabit/sec bit of 
communication channel through the existing infrastructure available with the 
optical community. Once such a quantum computer is available at remote site, 
these packets acting as "quantum software" can be transferred through high- 
speed communication channels. Thus, it is possible to carry out quantum 
computation at a remote distance with the proposed scheme of shaped pulses 
for terabit/sec communication and molecular control. For examples, see, 
Goswami D, Quantum Physics 0312192 (Dec. 24 2003); Garcia-Ripoll JJ, et al., 
Phys. Rev. Let., 91(15) 157901 (2003); Yang W, et al., J. Opt. Comm. Vol. 22, 
No. 1, pp. 694-697 (2001); Leibfried D, et al., J. Mod. Optic. 50 (6-7): 1115-1129 
(Apr-May 2003); Pazy E., et al,. Phys. Rep. 374 (6): 385-481 (Feb. 2003); 
Goswami D., Phys. Rep. 374 (6): 385-481 (Feb. 2003); Brattke S., et al., Phys. 
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Rev. Lett. 86 (16): 3534-3537 (Apr. 16 2001); and Weinacht T.C., et al, Nature 
397 (6716): 233-235 (Jan. 21 1999); in addition U.S. Patent No. 6,678,450 which 
issued to Franson on January 13, 2004; and US Patent 5,793,091 to Devoe 
issued August 11, 1998; all the above patents are incorporated by reference 
herein. 

Optical Coherence Tomography 
[0111] A preferred embodiment of the present invention uses a laser 
system 221 for laser excitation or ionization with Optical Coherence Tomography 
("OCT"). In general, Figure 13 illustrates the OCT application of system 221 
wherein there is a femtosecond laser 223, a laser beam shaper 229, a human or 
animal tissue specimen 241 , an optical gate 251 and an image 253. Laser 223 
emits a laser beam pulse shorter than 1 picosecond. Shaper 229 is made of 
three parts; two dispersive elements 255 which sandwich a phase mask element 
257. Shaper 229 essentially prevents multiphoton excitation which can damage 
the person's or animal's DNA, as will be discussed in more detail as follows. An 
unshaped laser beam pulse is used to gate the ballistic photons to render the 
image for tomography use. Optical gating can be accomplished by up- 
conversion in a frequency doubling crystal or with a kerr-gate in liquid carbon 
disulphide. The construction of system 221 as illustrated supposes transmission 
imaging; the same end result can alternately be accomplished with back 
scattered imaging. Image 253 could be viewed like an x-ray-type image of the 
internal organs of the human or animal specimen but without harmful three 
photon exposure. The use of the shaped pulse in OCT provides for an increase 
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in laser intensity for better imaging while preventing the damaging effects caused 
by multiphoton excitation of healthy tissue. The MIIPS and BPS processes 
discussed hereinafter can be advantageously used to activate different dyes and 
other compounds within a human or animal tissue, to achieve compound specific 
or functional OCT or microscopy. The pulse shaper is used to prevent three- 
photon and higher order nonlinear optical processes such as continuum 
generation. Higher order processes usually lead to sample degradation, and in 
the case of living samples to DNA damage. Suppression of three-photon 
transitions of four orders of magnitude has been achieved using the Mil and BPS 
methods and this suppression can be coupled with optimization of two-photon 
signal from living specimens. Alternatively, a fluorescent contrast agent can be 
administered so that pulses shaped using BPS selectively excite the fluorescent 
agent targeted towards malignant tumors. The embodiment is expected to 
achieve functional deep tissue imaging. 

[0112] Referring now to Figure 15, a system setup for functional 
imaging using BPS is shown. The tissue has been injected a fluorescent 
contrast agent that is preferentially absorbed by tumors. The fluorescent is a pH 
sensitive dye or derivatized quantum dots. No time grating is required. 

Photodvnamic Therapy 
[0113] Another embodiment of the present invention uses a system 
also shown as 221 for laser excitation or ionization with photodynamic therapy 
("PDT"). PDT is a treatment that involves the combination of visible light and a 
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photosensitizer. Each factor is harmless by itself, but when combined with 
oxygen, can produce lethal cytotoxic agents that can inactivate tumor cells. This 
enables greater selectivity towards diseased tissue as only those cells that are 
simultaneously exposed to the photosensitizer, light and oxygen are exposed to 
the cytotoxic effect. The dual selectivity of PDT is produced by both a 
preferential uptake of the photosensitizer by the diseased tissue and the ability to 
confine activation of the photosensitizer to this diseased tissue by restricting the 
illumination to that specific region. Therefore, PDT allows for the selective 
destruction of tumors while leaving normal tissue intact. 

[0114] PDT is based on the concept that (1) certain photosensitizers 
can be localized (somewhat preferentially) in neoplastic tissue, and (2) 
subsequently, these photosensitizers can be activated with the appropriate 
wavelength (energy) of light to generate active molecular species, such as free 
radicals and singlet oxygen ( 1 0 2 ) that are toxic to cells and tissues. PDT is a 
binary therapy, and a potential advantage of PDT is its inherent dual selectivity. 
First, selectivity is achieved by an increased concentration of the photosensitizer 
in target tissue, and second, the irradiation can be limited to a specified volume. 
Provided that the photosensitizer is nontoxic, only the irradiated areas will be 
affected, even if the photosensitizer does bind to normal tissues. Selectivity can 
be further enhanced by binding photosensitizers to molecular delivery systems 
that have high affinity for target tissue. For photoactivation, the wavelength of 
light is matched to the electronic absorption spectrum of the photosensitizer so 
that photons are absorbed by the photosensitizer and the desired photochemistry 



51 



6550-000057/CPE 



can occur. Except in special situations, where the lesions being treated are very 
superficial, the range of activating light is typically between 600 and 900 nm. 
This is because endogenous molecules, in particular hemoglobin, strongly 
absorb light below 600 nm and therefore capture most of the incoming photons. 
The net effect would be the impairment of penetration of the activating light 
through the tissue. The reason for the 900 nm upper limit is that energetics 
beyond this wavelength are insufficient to produce 1 0 2 , the activated state of 
oxygen, perhaps critical for successful PDT. 

[01 15] The photochemical and photophysical principles of PDT are well 
known in the art. Briefly, upon illumination, the photosensitizer is excited from the 
ground state (S 0 ) to the first excited single state (S^, followed by conversion to 
the triplet state (T^ via intersystem crossing. The longer lifetime of the triplet 
state enables the interaction of the excited photosensitizer with the surrounding 
molecules, and it is generally accepted that the generation of the cytotoxic 
species produced during PDT occurs whilst in this state. 

[01 16] The excited triplet state can react in two ways, defined as Type I 
and Type II mechanisms. A Type I mechanism involves hydrogen-atom 
abstraction or electron-transfer reactions between the excited state of the 
sensitizer and a substrate that is either biological, a solvent or another sensitizer, 
to yield free radicals and radical ions. These free radical species are generally 
highly reactive and can readily interact with molecular oxygen to either generate 
reactive oxygen species such as superoxide anions or hydroxyl radicals or can 
cause irreparable biological damage. These reactions produce oxidative damage 
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that is eventually expressed as biological lesions. By contrast, a Type II 
mechanism results from an energy transfer between the excited triplet state of 
the sensitizer and the ground-state molecular oxygen, generating the first excited 
state of oxygen, singlet oxygen. This zwitterionic species is extremely reactive 
and can interact with a large number of biological substrates, inducing oxidative 
damage and ultimately cell death. While it is generally accepted that Type II 
processes predominate during PDT and that singlet oxygen is the primary 
cytotoxic agent responsible for the biological effects displayed, Type I reactions 
become more important at low oxygen concentrations or in more polar 
environments. However, the initial reaction is of lesser importance as both Type 
I and Type II reactions lead to similar oxidative damage and comparable free 
radical chain-reactions in the presence of oxygen. The overall effect of either a 
Type I or Type II reaction pathway is the production of oxidative damage within 
the target cell that will ultimately lead to tumor destruction. 

[0117] Under special circumstances (short pulse, high intensities of 
irradiation), the upper excited states may be populated, and complex 
photophysical and photochemical processes may originate from these states, 
resulting in increased or decreased phototoxicity, which may include oxygen- 
independent mechanisms such as DNA mutation. 

[01 1 8] Photosensitizers are compounds that are capable of absorbing 
light of a specific wavelength and transforming it into useful energy. In the case 
of PDT, this would involve the production of lethal cytotoxic agents. There are 
hundreds of natural and synthetic dyes that can function as photosensitizers for 
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PDT, ranging from plant abstracts to complex synthetic macrocycles. The key 
characteristic of any photosensitizer is its ability to preferentially accumulate in 
diseased tissue and to then generate cytotoxic agents to induce the desired 
biological effect. (For examples of such photosensitizers See W.M. Sharman, et 
al., "Photodynamic therapy: basic principles and clinical applications," Drug 
Discovery Today 4(1 1):508-517 (1999); T. Hasan, et al., "Photodynamic Therapy 
Of Cancer," Chapter 40 in Holland Frei Cancer Medicine, BC Dekker Inc. (2003); 
W.M. Sharman, et al., "Targeted photodynamic therapy via receptor mediated 
delivery systems," Adv. Drug Delivery Rev. 56(1):53-76 (January 2004); and Roy 
I., et al., "Ceramic-based nanoparticles entrapping water-soluble photosensitizing 
drugs: A novel drug carrier system for photodynamic therapy." J. Am. Chem. 
Soc. 125:7860-7865 (2003)). 

[0119] In general, Figure 13 also illustrates the PDT application of 
system 221 , but optical gate 251 and image 253 are not required. Shaper 229 
allows two-photon excitations but essentially prevents three-photon excitation. 
Shaper 229 enhances the laser-induced activity of a therapeutic agent which 
prevents damage of healthy tissue. Use of laser beam pulse shaping of the 
present invention should provide superior control and results for PDT 
applications as compared to those practically possible with conventional methods 
as disclosed, for example, in U.S. Patent No. 6,042,603 entitled "Method for 
Improved Selectivity in Photo-Activation of Molecular Agents" which issued to 
Fisher et al. on March 28, 2000, and is incorporated by reference herein. 
Alternately, the pulse shaper can be tuned to target cancerous cells for 
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multiphoton gene therapy or destruction, with or without the presence of a 
therapeutic agent, without damaging healthy tissue. The Ml I PS and BPS 
processes discussed hereinafter can be advantageously used to activate only 
certain pharmaceuticals or chemicals, or used to allow the laser pulse to enter 
human or animal tissue to a known depth, based on the phase tuning and 
associated nonlinear spectrum tuning of the laser beam pulse. The pulse shaper 
is used to prevent three-photon and higher order nonlinear optical processes 
such as continuum generation. Higher order processes usually lead to sample 
degradation, and in the case of living samples to DNA damage. Suppression of 
three-photon transitions of four orders of magnitude has been achieved using the 
MM and BPS methods and this suppression can be coupled with optimization of 
two-photon signal from living specimens. 

[0120] General applications of lasers for biomedical purposes are well 
known for diagnostic tools, surgical tools and for imaging purposes. The 
immediate extension of ultrashort pulse shaping technology for biomedical 
applications holds many promises, to further necessitate looking into its domain 
of influence. Though very few applications for ultrafast pulse shaping technology 
are currently in use, there are very strong indications as to where it would lead. 
One of the most commonly adopted methods, for imaging, in recent times for 
three-dimensional profile measurement is optical coherence tomograph or a 
white light interferometer, which uses a broadband, low coherence light source. 
Recently the principle of femtosecond pulse shaping by spectral modulation has 
been used in conjunction with the joint transform correlator to make a spatio- 
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temporal joint transform correlator. The advantage of such a technique has been 
that it essentially removed the need of 1-D depth scanning and thereby avoided 
the long measurement times involved. Consequently, this eliminates the 
electronic computation needed to obtain the object image, and so it can be 
implemented as an all-optical set-up. Initially, this was demonstrated as a surface 
measurement set-up, however, as a natural extension, it was easily extended for 
providing tomographic sectioning of biological samples. In fact, with the use of 
principles of pulse shaping, a depth resolution of 70 ^m was achieved. 
Furthermore since there is no contact between the probe and the tissue, it is a 
useful non-invasive technique, which provides the physician with near- 
histological resolution imaging of sub-surface tissue morphology, potentially 
aiding in biopsy site selection and thus approaching the goal of "optical biopsy". 

Nonlinear Optical Excitation Spectroscopy 
[0121] In one embodiment of the invention, the goal is to develop laser 
systems capable of generating ultrashort pulses with unprecedented control over 
the spectral profile and phase properties of the pulses. Specifically, the systems 
are capable of generating ultrashort pulses that are within 1% of the transform 
limit, as determined by accurate pulse characterization, by a novel method that 
incorporates phase characterization and compensation in a single pulse 
shaping/characterization unit. The same unit will provide calibrated synthesis of 
arbitrary pulse shapes. One embodiment of the invention, a laser system forms 
an integral part of a microscope capable of functional imaging, a method 



56 



6550-000057/CPE 



whereby pulses are tailored to excite chromophores sensitive to their microscopic 
chemical environment (pH, Ca**, Na + gradients). 

[0122] In another embodiment of the invention, a laser system is 
amplified and shaped by a two-dimensional spatial light modulator. The system is 
capable of providing single-shot multiphoton excitation spectra over a 20-30 nm 
range, with 0.05 nm resolution. Two-photon excitation spectra of molecules are 
usually acquired point by point and are prone to order-of magnitude errors. The 
system is capable of acquiring 1000 spectra per second over 30 nm excitation 
regions, and referencing them to a standard. The system improves the accuracy 
and speed of nonlinear frequency resolved excitation spectroscopy and cross 
section measurements by orders of magnitude. This information is central to 
selection and utilization of nonlinear optical materials. 

[0123] Nonlinear optical excitation spectroscopy (NOES) involves the 
measurement of nonlinear processes such as two- or three-photon excitation 
cross-sections and nonlinear optical susceptibilities as a function of wavelength. 
These measurements are of great importance for the characterization of 
nonlinear optical materials and nonlinear optical chromophores such as laser 
dyes and quantum dots. Two-photon excitation spectroscopy, for example, 
requires a tunable laser source that is both narrow in frequency but highly 
efficient in multiphoton excitation. Typically, a narrow bandwidth tunable laser is 
used to obtain data from 700 to 950 nm, acquiring the data point by point. 
Nonlinear spectroscopy is extremely sensitive to transverse mode quality, 
spectral phase, characteristics of the focal spot, pulse-1to-pulse intensity 
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variations, and wave front deformations; all of which can cause higher or lower 
order processes that contaminate the result. All of these factors make point-by- 
point acquisition prone to systematic errors. 

[0124] A practical but naive approach would be to generate a broad 
bandwidth pulse and then to scan a slit across its spectrum to narrow the 
spectral resolution at specific wavelengths. This approach reduces the efficiency 
of second and higher order processes by two or more orders of magnitude. 

[0125] To solve these problems, an embodiment of the invention is 
based on binary phase masks based on the principles of intrapulse interference 
(MIIPS) and then optimized using learning calculations afforded by an 
evolutionary learning calculation. The binary phase mask and the corresponding 
output for transform limited and for shaped pulses are shown in Fig. 8, together 
with a demonstration of how these masks are used to tune the location where 
second order processes take place. The red trace is a simulation of the data. The 
laser is focused on a thin BBO crystal and the frequency-doubled output is 
collected with a spectrometer. Each phase mask is designed to produce a 
narrow and high contrast maximum at a different wavelength. These results were 
obtained with 18 fs pulses using a 128-pixel pulse shaper. Notice that the 
narrowed output is about 65% of the efficiency observed with TL pulses. Using 
an amplitude mask to block the spectrum of the laser, as shown in Fig. 11, 
above, resulted in only 3% conversion efficiency. 

[0126] An embodiment of the invention acquires spectra across a large 
bandwidth in a single laser shot. The high data throughput permits very accurate 
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calibration to some of the best-known standards such as the nonlinear optical 
crystal KDP and the laser dye Rhodamine 6G. The laser system is not be 
scanned; in fact, the amplitude of the laser field remains unchanged throughout 
the experiment. This permits outstanding calibration of all the factors that could 
affect nonlinear output. 

[0127] A diffractive two-dimensional programmable phase modulator 
originally designed by Hamamatsu for optical computing, is ideally suited for this 
task. The PPM X8267 is a 1024x768 pixel electrically-addressed phase 
modulator using an image transmitting element to couple and optically-addressed 
PAL-SLM (parallel aligned nematic liquid crystal spatial light modulator). The 
number of pixels available for pulse shaping is critical in the design. The 
generation of well-defined pulses requires the introduction of a well-defined 
phase functions. Pixilation effects, where a smooth phase is replaced by a 
jagged, under-sampled phase can be detrimental. The PPM units are illuminated 
by a spatial image, which is projected into the PAL-SLM. The projection system 
provides high accuracy reproduction without pixel borders, making it ideal for 
pulse shaping purposes. The femtosecond laser that is modulated is spectrally 
dispersed by a 300 lines/mm gold coated diffraction grating and collimated by a 
gold-coated cylindrical mirror forming a reflective 2f arrangement. The PAL-SLM 
is placed at the Fourier plane (where best spectral resolution is achieved). The 
reflected light is collected by an optical arrangement identical to the input one. 
Chromatic aberrations and group velocity dispersion are kept to a minimum by 
the all-reflective optical design. Nonlinear optical distortions and optical damage 



59 



6550-000057/CPE 

are minimized by the use of cylindrical rather than spherical optics. The 
spectrums are dispersed on the horizontal dimension. The pixels in the vertical 
dimension contain different phase masks that scan the wavelength where 
nonlinear excitation takes place. A schematic of the phase mask is shown in Fig. 
14. 

[0128] A embodiment of the invention, illustrated in Fig. 15, produces 
amplified sub-20 fs pulses (~10 6 times more energy per pulse than laser System 
A), but at a repetition rate of 1 kHz. The bandwidth of the pulses is 50 nm 
FWHM, and is used to obtain NOES with about 0.2 nm resolution. The seed 
pulses are produced by a K&M Labs oscillator which is pumped by a Verdi laser 
(Coherent). The amplifier is an Odin, multipass amplifier from Quantronix, 
pumped by their Quantronix Darwin Q-Switched Nd:YLF, capable of delivering 1 
mJ per pulse at 1 kHz. The Quantronix system is quoted for sub-35-fs 
performance. However, the amplifier is capable of preserving most of the 
oscillator bandwidth. The amplifier design is based on the multi-pass platform 
that Kapteyn and Murnane used to demonstrate generation of 17 fs pulses. To 
achieve the very short pulse duration a grating and prism compressors are used. 

[0129] The laser system is installed on a 5' x 12' vibrational^ isolated 
optical bench (Newport). Pulse characterization is carried out by FROG, SHG- 
FROG, and MIIPS. For FROG, an optical delay line (Aerotech) is used. Data 
collection is accomplished using an Ocean Optics spectrometer, a Boxcar 
Averager/Integrator (SRS), a 500 MHz digital oscilloscope (Infinium, Hewlett 
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Packard), and controlled with a personal computer running LabView, with a GPIB 
IEEE 488 controller. 

[0130] Two-dimensional pulse shaping allows single-shot NOES data 
to be acquired from nonlinear optical materials, two-photon excitation spectra of 
laser dyes, biologically tagged fluorescent markers, quantum dots, optical 
switches and other optical materials. The wavelength range is limited to the 700 
to 900 nm window by the pump pulses. The laser system is used to pump a non- 
collinear optical parametric amplifier, producing tunable broad bandwidth pulses 
from 450 to 1 300 nm. These pulses are directed to the two-dimensional pulse 
shaper when characterization of a material is required at very different 
wavelengths. 

[0131] The MM PS method, described previously, requires the scanning 
of a reference phase across the spectrum of the laser. This usually requires 
acquisition of the SHG output spectrum for -100 different phase function 
positions. To achieve this task in a single shot, the vertical dimension is divided 
into 128 different sections. Each section, containing 6x1024 pixels, contains the 
reference phase function. The key is that the phase function is displaced in each 
of the sections to provide the data as a function of 6 required for MIIPS. Once the 
femtosecond laser pulse undergoes the two-dimensional phase modulation it is 
focused with a short focal length -100 mm cylindrical mirror on a thick SHG 
crystal. On the horizontal axis, the thick SHG crystal causes a spectral 
dispersion that is analogous to the one used in the GRENOUILLE. On the 
vertical axis, one obtains different sets of sets of SHG spectra because of the 
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different phase functions encoded by the PAL-SLM. The entire two-dimensional 
image, containing the one-shot Ml IPS data is then imaged onto a CCD for 
collection and analysis. 

Communications 

[0132] With reference to Figure 16, another preferred embodiment of a 
laser excitation system 421 of the present invention employs a femtosecond 
laser 423, an optical fiber 451 , a laser beam pulse shaper device 429, a laser 
beam pulse un-shaper device 453, and a receiver 441 which includes an optical 
switch or sensor and the related circuitry and electrical control unit. Laser 423 
emits a series of laser beam pulses, each shorter than 1 ps, into the connected 
fiber 451 . Pulse shaper device 429 is of a predetermined mask type with a fixed 
pulse characteristic varying shape (such as with calculated sine wave surface 
shapes) and has three elements connected to fiber 451: a dispersive element 
455 such as a fiber that incorporates a diffraction grating; a phase mask element 
457 that can be made using a doped glass or polymer sheet; and a dispersive 
element 459, like element 455 but reversed, for accepting spectrally dispersed 
light and coupling it back to fiber 451 . 

[0133] The shaped laser beam pulse is capable of traveling long 
distances through fiber 451 without suffering nonlinear distortion because of the 
unique phase function imprinted or formed on shaper device 429. For example, 
the red color spectrum may be advanced in front of the blue color spectrum in a 
precise sine manner. Un-shaper device 453 subsequently reverses the phase 
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changes introduced by shaper device 429. It is constructed the same as the 
shaper device but with a different phase mask element 461 that compensates for 
the pulse characteristic changes made by mask element 457. Alternately, an 
acousto-optic modulator or transient grating can be used for optical switching 
through constructive or destructive reference of waves. Shaping and unshaping 
can also be accomplished by means of a chirped mirror or spectral masks. 

[0134] Thus, the present invention's ability to precisely control the laser 
beam pulse shape or other characteristic, especially for nonlinear or multiphoton 
emissions, significantly improves the quality of the communication transmission 
while minimizing self-focusing, self phase modulation and possible destruction of 
the fiber. The pulse characteristic control of ultrafast laser beam pulses, as 
described in all of the embodiments herein, should minimize, if not prevent, 
multiplicative noise effect disruption of nonlinear propagation channels in fiber 
optic lines, as discussed in Mitra, et al., "Nonlinear Limits to the Information 
Capacity of Optical Fibre Communications," Nature, vol. 411, pp. 1027-1030 
(June 28, 2001). It is further envisioned that this type of pulse shaping system 
can be employed within salt water oceans for submarine-to-submarine 
communications using short laser pulses of 1 ps or less. This type of pulse 
shaping can be used to induce solution formation to achieve minimally distorting 
pulses for communications. Moreover, MIIPS can be used to measure the 
distance of a fs laser emitter by determining the magnitude of the acquired 
second order phase modulation as the laser pulse transmits through air or water. 
This method does not require echo or reflection. In water longer pulses (1 ps) 
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are desired because of the much greater dispersion. Depending on the 
transmission medium, air or water, and the distances expected different pulses 
are required. For air, short pulses with durations between 10-20 fs will be 
preferred. For water, pulses with much longer durations will be preferred, for 
example for 100 m distance 100 ps pulses would be preferred. 

[0135] Using the BPS method, the data density that can be achieved 
per pulse is equal to or less than the number of pixels in the SLM divided by 2. 
For example, in an embodiment in the SLM 429 has 256 pixel resolution, the 
maximum data density is 256/2 or 128 bits per pulse. In this embodiment, if the 
pulse rate is 100 MHz, the data rate or bandwidth would be 10 10 bits per second. 
Other embodiments of the invention vary the data rate or bandwidth based on 
pixel resolution and/or laser pulse rate. As laser pulse rates increase and as 
SLM pixel resolution increases, the bandwidth achieveable by this invention will 
increase. It is advantageous that only a nonlinear optical detector would be able 
to decode the signal and, in contrast, a simple light detector would not be able to 
determine any information from the pulses. The fact that it can be used 
asynchronously is ideal for mobile or distant, and intermittent communications. 

[0136] Referring to Figures 17, 18A and 18B, a fourth preferred 
embodiment of the system of the present invention is used for fiber optic 
communications. Multiple transmission users who are each sending a 
communications message or signal are using a communications device such as 
a telephone 491, personal computer, facsimile machine or the like, at remote 
locations from each other. These remote transmitters are connected to a "smart" 
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main transmitter assembly which includes a computerized, central processing 
unit 493 through electric wires, fiber optic cables, microwave signals or the like. 
A phase modulated pulse shaper 505 is actively controlled by CPU 493. Laser 
509 and shaper 505 are also contained as part of the main transmitter assembly. 
Laser 509 emits an ultrashort laser pulse which is carried within a fiber optic 
cable 497 after shaping. The ultrashort laser beam pulses have a duration of 
about 100 femtoseconds based on currently available fiber optic cable limitations 
but pulse durations of less than 50 femtoseconds would be preferred and those 
of 10 or less femtoseconds would be the most desired if fiber optics allow for 
such in the future. For example, photonic band gap materials such as optical 
fibers with holes therein may allow for use of approximately 10 femtosecond 
pulses. 

[0137] Pulse shaper/phase mask 505 encodes each laser beam pulse 
phase, using a binary phase mask. The second harmonics contains multiple 
peaks, by way of example, but not limitation, in the frequency domain, thus 
revealing the encoded message. This allows encoding of routing addresses and 
the associated communications information to be encoded within each laser 
beam pulse based on CPU control of the laser beam emissions in combination 
with actively varied shaping of each emitted pulse. 

[0138] A "dumb" central receiver 501, one that does not require an 
additional laser or complex computational capabilities, is connected to the 
downstream end of fiber optic cable 497. Receiver 501 includes a focusing lens 
515, a thick SHG crystal 507' and a detector 511. Each laser beam pulse 
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transmitted through fiber optic cable 497 is dispersed onto lens 515 which serves 
to focus and direct each pulse, in a converging angular manner, onto crystal 
507'. A thick optical crystal 507' is defined herein as one having a transmissive 
path thickness of greater than about 0.5 millimeters while a thin optical crystal 
507 (see Figure 15) is defined herein as having a transmissive path thickness 
less than about 0.5 millimeters. The preferred thickness for the thick crystal is 
approximately 3.0 millimeters for 50 femtosecond or less pulse duration and 5.0 
millimeters for a 50 to 200 femtosecond pulse duration. Thick crystal 507' 
creates a second order harmonic and second order spectrum within each pulse 
as previously shaped by the pulse shaper. In other words, the thick crystal 
disperses essentially the entire color spectrum without use of a separate 
spectrometer because of the phase matching angle requirement. 

[0139] Each separated color frequency angularly dispersed from the 
thick crystal is encoded by the pulse shaper to contain individual communication 
routing addresses and the actual communications information, which is then 
detected by a multiplexer-type of detector 51 1 such as a CCD camera employing 
a linear array. Alternately, detector 511 is a two-dimensional array that can be 
used to achieve higher data densities by adding one more dimension. It is also 
alternately envisioned that detector 511 is an array of optical fibers that are 
connected to remote controllers/sub-detectors. The data can be read 
asynchronously using only the transmission pulse containing the information and 
not additional reference pulse. A single detector 51 1 is operable to digitize the 
detected signals carried in each pulse as separated through the spectrum and 
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transmit them through wires, fiberoptics, microwaves or the like to individual 
decoding microprocessor controllers 503 within or external to receiver 501 . A set 
of prestored variables or dencryption information or key is located within memory 
of each controller 503 in order to decode each corresponding digitized 
communication signal received by detector 511 without requiring synchronous 
communication transmissions (in other words, a second laser pulse that provides 
a complimentary phase) from transmitter 495. The decoded communications are 
then sent to the end users who receive such by telephones 505, personal 
computers, facsimile machines or the like at the identified routing addresses 
desired. Alternately, controllers 503 can be replaced by simple light detection 
devices such as photodiodes which can be employed in a digitized on/off self- 
switching mode based on the signal detected by detector 51 1 to control or send 
information to remote destinations. It is significant that interferometry and 
synchronous laser pulses are not required for decoding the transmitted 
information with the presently preferred communications embodiment of the 
present invention. It is also noteworthy that pulse shaper 505 can encode each 
pulse by use of second harmonic generation or any other non-linear mixing 
method including, but not being limited to, frequency mixing, difference frequency 
mixing, and four wave mixing. 

[0140] The present invention should be contrasted to a prior 
experiment which employed a difficult and a synchronous reference pulse at the 
decoder for supplying a complimentary phase to control the emission of a single 
specific wavelength. This is disclosed in Z. Zheng and A. Weiner, "Coherent 



67 



6550-000057/CPE 

Control of Second Harmonic Generation Using Spectrally Phase Coded 
Femtosecond Waveforms," Chemical Physics 267, p. 161 (2001); this prior 
approach, however, required pulses which overlap in time and space, which is 
difficult to control, and only for a single pulse frequency. 

[0141] While various embodiments have been disclosed herein, it 
should be appreciated that other modifications may be made that are covered by 
the system and methods of the present invention. For example, alternate lasers, 
chemicals, optics and computer controllers can be employed as long as they 
function as described. The description of the invention is merely exemplary in 
nature and, thus, variations that do not depart from the gist of the invention are 
intended to be within the scope of the invention. Such variations are not to be 
regarded as a departure from the spirit and scope of the invention. 
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